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Abstract 

Measurements of inclusive jet, dijet and three-jet cross sections in photon-proton 
interactions are presented. These measurements provide new tests of QCD, 
constrain the parton densities of the photon, and allow searches for new physics. 
Measurements of the mean subjet multiplicity are reported and used to test the 
differences between quark and gluon jets. 



1. Introduction 

Measurements of cross sections for processes 
involving a large momentum transfer (e.g. jets 
in 7p interactions) are compared to next-to- 
leading order (NLO) calculations to test QCD. 
These measurements also provide a test of the 
parametrisations of the photon parton densities and 
may be used in global analyses to constrain the 
parton distributions. 

The investigation of the internal structure of jets 
gives insight into the transition between a parton 
produced in a hard process and the experimentally 
observable spray of hadrons. Measurements of jet 
substructure allow the study of the characteristics 
of quark and gluon jets. 

2. Jet cross sections in interactions 

At HERA, positrons of energy = 27.5 GeV 
collide with protons of energy Ep — 820 GeV. The 
main source of jets at HERA is hard scattering 
in 7p interactions in which a quasi-real photon 
{Q^ ~ 0, where is the virtuality of the photon) 
emitted by the positron beam interacts with a 
parton from the proton to produce two jets in the 
final state. At leading order (LO) QCD, there are 
two processes which contribute to the jet production 
cross section: the resolved process (figure |l|a) in 
which the photon interacts through its partonic 
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content, and the direct process (figure Qb) in which 
the photon interacts as a point-like particle. 

The cross section for jet production at LO QCD 
in 7p interactions is given by 
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where f^/eiv) is the flux of photons in the 
positron, usually estimated by the Weizsacker- 
Williams approximation (y is the fraction of the 
positron energy taken by the photon); fj/pixp^fij^) 
are the parton densities in the proton, determined 
from e.g. global fits ||^ (xp is the fraction of the 
proton momentum taken by parton j and /ii? is the 
factorisation scale); and d<7('j(i)j —^ jet jet) is the 
subprocess cross section, calculable in perturbative 
QCD. In the case of resolved processes, there is an 
additional ingredient: fi/'yix.y, fijr) are the parton 
densities in the photon, for which there is only 
partial information (x-y is the fraction of the photon 
momentum taken by parton i). The integrals are 
performed over the phase space represented by 

2.1. The iterative cone algorithm 

In hadronic type interactions, jets are usually recon- 
structed by a cone algorithm Experimentally, 
jets are found in the pseudorapidity (77) — azimuth 
{ip) plane using the transverse energy flow of the 




event. The jet variables are defined according to 
the Snowmass Convention Q, 

„jet _ Y^jET-Vi . ,iet _ IZjET-Vi 

In the iterative cone algorithm, jets are found by 
maximising the summed transverse energy within a 
cone of radius R. 

3. Inclusive jet cross sections 

Inclusive jet cross sections have been measured |^ 
using the 1995-1997 ZEUS § data (which amounts 
to an integrated luminosity of £ ^ 43 pb~^) as a 
function of the jet transverse energy. The jets have 
been searched for with an iterative cone algorithm 
with R = 1. The measurements have been 
performed for jets of hadrons with Ei^^ between 
17 and 74 GeV and 77^*=* between -0.75 and 2.5, 
and are given for the kinematic region defined by 
0.2 < y < 0.85 and < 4 GeV^. 

Figure ||a shows the measured da/dEi^^ (black 
dots). The systematic uncertainties not associated 
with the absolute energy scale of the jets have 
been added in quadrature to the statistical errors 
(thick error bars) and are shown as thin error bars. 
The shaded band represents the uncertainty on the 
energy scale of the jets. The data show a steep fall- 
off over four orders of magnitude in the measured 
range. 

3.1. NLO QCD calculations 



There are several complete calculations of jet cross jets |16 
sections at NLO for interactions |^,^,^,|l^]. Two 



mainly in the treatment of the real corrections. A 
detailed comparison of the calculations |11| show 
that for NLO calculations of dijet cross sections 
and LO calculations of three-jet cross sections a 
reasonable agreement is found, and the differences 
are up to ^ 5%. 

The curves in figure ^a are NLO QCD 
calculations [Q, ^ using different parametrisations 
of the photon structure function: GS96 [Q (solid 
line), GRV-HO P] (dashed line) and AFG @ 
(dot-dashed \iney The CTEQ4M [|| proton 
parton densities have been used in all cases. In 
the calculations shown here, the renormalisation 
and factorisation scales have been chosen equal 
to Ei^ and as was calculated at 2-loops with 
= 296 MeV. The NLO calculations give a 
reasonable description of the data. Figure ||b shows 
the fractional differences between the measured 
da/dEi^* and the NLO calculations based on GRV- 
HO. 



3.2. Comparing theory and experiment 

To perform tests of QCD and to extract information 
on the photon parton densities, the experimental 
and theoretical uncertainties must be reduced as 
much as possible. 

Among the main experimental uncertainties the 
presence of a possible underlying event, which is 
the result of soft interactions between the partons 
in the photon and proton remnants (for resolved 
events), and is not included in the calculations. 
The uncertainty of the measurements due to the 
underlying event is reduced by decreasing the cone 
radius or by increasing the transverse energy of the 



types of corrections contribute at NLO: the virtual 
corrections which include internal particle loops and 
the real corrections which include a third parton 
in the final state. The existing calculations differ 



On the theoretical side, since calculations are 
made only at NLO, the implementation of the 
iterative cone jet algorithm in the theory does not 
match the experimental procedure exactly. The 
theoretical uncertainty coming from this effect is 
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Figure 2. Inclusive jet cross section as a function of E-'r^ using the iterative cone algorithm. NLO QCD 
calculations are shown for comparison. 



reduced by using the longitudinally invariant fej- 
cluster algorithm |17|]. 



3.3. The kx cluster algorithm 

In the inclusive fcy cluster algorithm jets are 
identified by successively combining nearby pairs 
of particles until a jet is complete. The kx 
algorithm allows a transparent translation of the 
theoretical jet algorithm to the experimental set-up 
by avoiding the ambiguities related to the merging 
and overlapping of jets and it is infrared safe to all 
orders. 

Figure ^ shows da/dE^* for jets found using 
the kr cluster algorithm. The NLO QCD 
calculations, using the current knowledge of the 
photon structure, are able to describe the data 
within the present experimental and theoretical 
uncertainties. 

A comparison between da/dE^* measured using 
the cone and the kr algorithms has been made (see 
figure ^). The differences between the measured 
cross sections are typically smaller than 10%. This 
comparison shows that the cone with R = 1 and 
the fc-r cluster algorithms probe the underlying 
parton dynamics in a comparable way. Therefore, 
in the experiment, the choice of jet algorithm is 
not crucial. The use of the fc^ cluster algorithm 
is dictated by the need to reduce the theoretical 
uncertainties. 



4. Jet substructure 

The use of the fc^ cluster algorithm allows the study 
of the internal structure of the jets in terms of 
subjets. Subjets are jet-like objects within a jet and 
are resolved by reapplying the kr cluster algorithm 
until for every pair of particles i and j, 



dij = min{ETt, ETj)'^[ir]i - Vi)^ + {^i - 



> VcntiE^T 

Measurements of the mean subjet multiplicity (< 
«subjct >) have been performed by ZEUS ||l8| using 
an inclusive sample of jets with Ei^j^* > 15 GeV and 
— 1 < 77^*^* < 2 in the kinematic range defined by 
0.2 < y < 0.85 and < 1 GcV^. 

Figure ^ shows < ngubjot > as a function 
of 2/cut. The data (black dots, the statistical 
and systematic uncertainties are included but are 
smaller than the dots) show that < ngubjot > grows 
as ?/cut decreases within the measured range. The 
lines are calculations from the leading-logarithm 
parton-shower Monte Carlo's PYTHIA 119] and 
HERWIG |§. The calculations based on PYTHIA 
give a good description of the data. Figure ||b 
shows < risubjet > as a function of rj^^* for ycut = 
0.01: < risubjot > increases as ry-''^* increases. The 
comparison with the predictions for quark and 
gluon jets shows that the increase in < risubjct > 
as ry-''^' increases is consistent with the predicted 
increase in the fraction of gluon jets. 
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Figure 3. Inclusive jet cross section as a function of using the kx cluster algorithm. NLO QCD calculations 
are shown for comparison. 
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Figure 4. Ratio between the measured da/dE!^ using the iterative cone and the kr cluster algorithms. 



5. Dijet cross sections 

Dijet cross sections have been measured |Q using 
the 1994 - 1997 HI ^ data (which amounts to 
an integrated luminosity of £ ~ 36 pb~^) as a 
function of the transverse energy of the leading 
jet and the average transverse energy of the two 
leading jets. The jets have been found using the 
kx cluster algorithm. The measurements have been 
performed for jets of hadrons with Ei!^'^^ > 25 GeV, 

S^'^*^ > 15 GcV and -0.5 < tj^"^ < 2.5, and are 
given for the kinematic region defined by y < 0.9 
and Q2 < 4 GeV2. 



Figure ^ shows the measured cross sections (black 
dots). The systematic uncertainties (including that 
associated with the absolute energy scale of the jets) 
have been added in quadrature to the statistical 
errors (inner error bars) and are shown as the 
outer error bars. The data show a steep fall-off of 
three orders of magnitude in the measured range. 
The histograms are the calculations using PYTHIA 
which provide a good description of the shape of the 
measured distributions. 
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Figure 5. Mean subjet multiplicity as a function of (a) ycut for —1 < ry-"^' < 2 and (b) 77-''^* for j/cut ~ 0.01. 
PYTHIA and HERWIG Monte Carlo calculations are shown for comparison. 
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Figure 6. Dijet cross sections as a function of (a) the transverse energy of the leading jet and (b) the average 
transverse energy of the two leading jets. PYTHIA Monte Carlo calculations are shown for comparison. 



6. High-mass dijet cross sections 

The dijet mass distribution M'''^ is sensitive to 
the presence of new particles or resonances that 
decay into two jets. The distribution of the angle 
between the jet-jet axis and the beam direction 
in the dijet centre-of-mass system (cosd*) reflects 
the underlying parton dynamics and is sensitive to 
the spin of the exchanged particle. New particles 
or resonances decaying into two jets may also be 
identified by deviations in the measured cos 9* 
distribution with respect to the QCD predictions. 

The cross section as a function of the dijet 
invariant mass has been measured by HI |2l| using 



the kx cluster algorithm between 48 and 148 GeV 
(figure |7|) . The data show a steep fall off of more 
than two orders of magnitude within the measured 
range. The histogram is the calculation using 
PYTHIA which gives a reasonable description of 
the shape of the measured distribution. There are 
large uncertainties in the normalisation of the LO 
QCD calculations which indicate the need for NLO 
corrections. 

High-mass dijet cross sections have been mea- 
sured H using the 1995 - 1997 ZEUS data as 
a function of M'^'^ and \cos9*\. The measure- 
ments have been performed for M'^'^ > 47 GeV and 
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Figure 7. Dijet cross section as a function of the dijet invariant mass. PYTHIA Monte Carlo calculations are 
shown for comparison. 
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Figure 8. High-mass dijet cross sections as a function of (a) the dijet invariant mass and (b) the dijet angular 
distribution. NLO QCD calculations are shown for comparison. 



cos6'*| < 0.8 using the fc-r cluster algorithm. 
The data show a steep fall-ofF in M'^'^ of 
three orders of magnitude in the measured range 
(figure ^). The measured da /d\ cos 9* \ rises as 
I cos 6** I increases (figure The NLO QCD 

calculations give a reasonable description of the 
measured distributions. The calculations based on 
GRV-HO are closer in magnitude to the measured 
cross sections. No significant deviation between 
data and NLO calculations is observed in the 
measured range of M'^-^ and | cos 6**1. 



7. Three-jet cross sections 

Measurements of three-jet cross sections provide a 
test of QCD beyond LO and allow a search for new 
phenomena. NLO calculations for three-jet cross 
sections in jp interactions are not yet available. The 
calculations shown here are LO for these processes 
and are therefore subject to large renormalisation 
and factorisation scale uncertainties. The cross 
section for three-jet production at LO is given by 



[ft/jix-^, da{-f\i]j jet jet jet). 



Figure 9. Three-body rest frame. 



Five parameters are necessary to uniquely 
determine the three-body phase-space. These are 
the three-jet invariant mass (M^j); the energy- 
sharing quantities and X4 (the jets are 
numbered 3, 4 and 5 in order of decreasing energy), 



X, 



_ 2Ei 



the cosine of the scattering angle of the highest 
energy jet with respect to the beam, 



\Pb\\P3\^ 

and "03, the angle between the plane containing 
the highest energy jet and the beam and the plane 
containing the three jets. The latter is defined by 



cos -03 



(P3 XpB)-(p4Xp5) 
|P3 XPB I |P4Xp5 I 



The definition of the angles 6*3 and V'3 is illustrated 
in figure 

Since 9^ involves only the highest energy jet, the 
distribution of cos^s in three-jet events is expected 
to follow closely the distribution of cos^* in dijet 
events. The ip3 angle, on the other hand, reflects 
the orientation of the lowest energy jet. 

Figure |lO| shows the three-jet cross section as 
a function of the transverse energy of the lowest 
energy jet and the three-jet invariant mass. The 
comparison to the calculations using PYTHIA 
Monte Carlo shows that the parton shower models 
describe the shape of the measured cross sections. 

The three-jet invariant mass cross section was 
measured by ZEUS in the kinematic region 
defined by | cos 6*31 < 0.8 and X3 < 0.95 (see 
figure |ll|). The curves in figure |ll| are the 0{a^) 
QCD calculations (Hl^ using the GRV-LO (l| 



parametrisations of t he p hoton structure function 
and the CTEQ4 LO llq| proton parton densities. 
The renormalisation aSd factorisation scales have 
been chosen equal to EJ^^^, the largest of the Ei!^* 



values of the three jets, as was calculated at 1- 
loop with = 181 MeV. The calculations give 

a good description of the data, even though they 
are LO for this process. Monte Carlo calculations 
are also compared to the data: they provide a good 
description of the data in shape, but the magnitude 
is 30 - 40% too low. 

The search for new particles or resonances 
decaying into two jets can be extended by looking 
for deviations in the distributions of the dijet 
invariant masses in three-jet events with respect 
to the predictions of QCD. Figure |l^ shows the 
dijet invariant mass distributions in three-jet events 
for all possible pairs of jets. The histograms are 
the predictions from the QCD-based Monte Carlo 
models. No significant deviation between data and 
calculations is observed up to the highest invariant 
mass value studied. 

Figure [l3| a and |l3| b show the X^ and 
distributions measured by ZEUS for M^j > 
50 GeV, I cos 6*3 1 < 0.8 and ^3 < 0.95. Calculations 
from different models are compared to the data: a 
pure phase space calculation does not describe the 
data. The 0{aal) QCD calculations are in good 
agreement with the data. 

The angular distribution of the lowest energy 
jet in three-jet events is a distinct probe of 
the dynamics beyond LO. The measured -03 
distribution (figure ^3|d) is drastically different from 
pure phase space and is in agreement with the 
0{aal) QCD calculations. The comparison of the 
data to parton shower models shows that the data 
favour color coherence. 

The measured cos ^3 distribution (figure [l^) 
indicates that the highest energy jet tends to 
go either forward (proton direction) or towards 
the rear (photon direction). The 0{aal) QCD 
calculations and the Monte Carlo models are in 
good agreement with the data. 
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Figure 10. Three-jet cross section as a function of (a) the transverse energy of the lowest-energy jet and (b) 
the three-jet invariant mass. PYTHIA Monte Carlo calculations are shown for comparison. 
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The variable x'^^^ , 

nBS ^ ^^.tsE^T^ 

gives the fraction of the photon energy invested in 
the production of the three-jet system and can be 
used to define resolved and direct processes in a 
meaningful way to all orders. 

Figure [ij shows the cos^a distribution for 
different regions of x^^^ . The data indicate that 
the highest energy jet tends to go in the forward 
direction for x^^^ < 0.5 (where LO resolved 
processes dominate) and it goes more backward as 



x^^^ increases. The 0{aal) QCD calculations 
are in good agreement with the data except for 
0.9 < x^^^ < 1. 

8. Conclusions 

Significant progress in comparing measurements 
and QCD calculations of jet cross sections in 7p 
interactions has been achieved; experimental and 
theoretical uncertainties have been reduced. 

NLO QCD calculations of inclusive jet and dijet 
cross sections and 0{aa'^) QCD calculations of 
three-jet cross sections describe reasonably well 
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Figure 12. Dijet invariant mass cross sections in three-jet events. PYTHIA Monte Carlo calculations are shown 
for comparison. 



the measurements. No significant deviation with 
respect to QCD predictions has been observed 
writhin the measured range of the variables studied. 
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